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Abstract - We reveal and study the topological transition in a metamaterial formed by par¬ 
allel nanowires of polaritonic material. When the dispersion transits from the elliptic (epsilon¬ 
positive) to hyperbolic (epsilon-indefinite) regime, a very specific isofrequency surface arises 
which implies an extraordinary Purcell factor in spite of noticeable optical losses. 

I. Introduction 

Wire media - an important class of hyperbolic metamaterials - were reviewed in d, where it was noticed 
that their key feature is strong spatial dispersion. The spatial dispersion distorts the isofrequency surfaces of 
regular wire media making them qualitatively different from hyperbolic ones - either flattened or locally deviating 
their general shape in the vicinity of the T point of reciprocal lattice. In the available literature, wire media 
are characterized by an indefinite complex permittivity tensor - a uniaxial tensor whose components have either 
positive or negative sign of the real part, whereas the imaginary part of these components keeps reasonably small. 
Of course, both axial and transverse components of this tensor may be also negative, but it is not interesting since 
the wire medium (WM) becomes opaque. More interesting is to notice that both these components can be positive 
and the WM may become an anisotropic dielectric with elliptic type of dispersion. 

Below we show that both hyperbolic and elliptic regimes are achievable with nano wires prepared of so-called 
polaritonic material, such as LiTa03. In such WM, the transverse component of the effective permittivity tensor 
keeps positive whereas the axial component (its real part) changes the sign over the frequency axis in the mid IR 
range. This means that there is a frequency at which the hyperbolic dispersion regime transits to the elliptic one and 
vice versa. In accordance to Il2|, this regime is called topological transition. Topological transition in a hyperbolic 
metamaterial implemented as a stack of metal-dielectric bilayers is already known O. On condition of low optical 
losses this regime is characterized by extraordinary density of photonic states El El- This density of states implies 
the so-called Purcell effect, and one of the most important implications of such a regime is a high Purcell factor - 
enhancement of the decay rate of a quantum emitter submerged into such a medium. Equivalently, this factor can 
be defined as the enhancement of the radiated power of a monochromatic subwavelength dipole source H . 

In this paper, we show that the radiation of a short dipole is strikingly enhanced in spite of rather high optical 
losses in our WM which would have obviously suppressed the Purcell effect in a usual hyperbolic medium. The 
reason of high Purcell factor is the combination of elliptic and the hyperbolic isofrequencies in the reciprocal space. 
The two branches are connected at two points where Re(5||) = 0. This becomes possible because the spatially 
dispersive effective permittivity depends on the axial component of the wave vector. As a result, neither at high 
nor at low spatial frequencies there is a cut-off: in the effective-medium approximation the spatial spectrum is 
unbounded. Another feature is the high Purcell factor for an axial dipole - in recent studies published after 01 the 
radiation of an axial dipole is suppressed, and the Purcell effect holds only for a transversely polarized dipole. 

II. Theory 

Optical properties of a WM sketched in Eig.l^in the effective-medium model is described by a uniaxial dyad 


e = 5 _lXoXo + ^uyoyo + ^||ZoZo, 


( 1 ) 
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Fig. 1: Long parallel wires (^r) which are set in a square lattice, b and a are the lattice constant and the radius of a 
dielectric rod, respectively, and 5h represents the host medium. 


where the transverse (5 _l) and axial (5||) components, in case of thin wires performed of negative-epsilon material, 
are given by the well-known spatially-dispersive model by Silveirinha (see e.g. in my. 


s± — 1 -|- 
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( 2 ) 


Here, Sr is the relative dielectric constant of the wires, /v denotes the volume fraction, ko is the free-space wave 
number and p represents the axial component of the wave vector. Also, is called the plasma wave number ex¬ 
pressed as /cp = (1/6) \/27r/[0.5275 + In (6/27ra)] (a is the radius of the wires and b denotes the lattice constant). 
Equation [^assumes that wires are embedded in free space. What is important in Eq.j^is that the axial component 
of the permittivity depends on the component of the wave vector parallel to the rod wires. In other words, the wire 
medium is spatially dispersive. 

In wire medium since it is a kind of uniaxial anisotropic medium, the transverse component of the wave vector 
is written as = s±kQ — which corresponds to ordinary or TE-waves, and = e\\[kQ — that 

corresponds to extraordinary or TM-waves. In this work, only TM-waves are important because in the vicinity 
of the short dipole, the reactive electromagnetic power stored as the evanescent waves are related to transverse 
magnetic (TM) waves. Based on the formula given above, if we assume that is larger than zero, TM-waves can 
propagate within the medium if either 1) k^ — f3‘^l£± > 0 and ey > 0, or 2) k^ — f3‘^l£± < 0 and < 0. Both 
conditions can be implemented at a single frequency: the first case resulting in an elliptic dispersion corresponds 
to ^ < koy^ and the second case giving the hyperbolic dispersion -to ^ > ko^^^. 

WM grants both these possibilities since its axial component of the effective permittivity (ey) depends on /3. 
Assume that for [3 = kQ^Je±. at a certain frequency, ey gets zero. After some algebraic manipulations based on 
Eq. i we can conclude that er = 1 + ^p/[/v(^o ~ ~ ^p)] ~ 1 ~ l//v < 0- Now, for P > ko^/^, it can be 

shown that ey < 0, for {3 < kQy/e]^, ey < 0, and at /3 = ey = 0. In other words, in wire medium we can 

join hyperbolic and elliptic dispersions together at the point where ^y = 0 and this is possible when the material 
of nanowires has negative permittivity er. 

To confirm our expectation on the Purcell factor, we have done several full-wave simulations using the CST 
software. For example, the fraction ratio fy = 0.0804 may correspond to a = 32 nm and b = 200 nm. Therefore, 
the value er ~ —11.4 is that allowing ey = 0 at a point p = This value for er corresponds to the 

permittivity of lithium tantalate (LiTa03) at the frequency / = 39 THz. Since the isofrequency contour in the 
plane {q — /3) is stretched along the axis (3, the high Purcell factor should be observed for a dipole creating 
the spatial spectrum which broadens over the optical axis than in the transverse plane. This implies the axial 
polarization of the dipole. Though in work 01 the high Purcell effect in WM is claimed for a transverse dipole 
producing more TM-polarized spatial harmonics, for our purposes, the isofrequency shape turns out to be more 
important, and the Purcell effect is more pronounced for an axial dipole. 
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Fig. 2: (a)- Dipole location and orientation in a finite-sized cubic sample of WM. (b)-Radiation enhancement of 
the dipole for three different sizes of the sample: 12x12, 16x16 and 20x20 nanowires with period b =200 nm. 
The blue dashed curve corresponds to the radiation enhancement of a transversal dipole in an infinite PEC (perfect 
electric conductor) wire medium with same geometry. 


We simulated an electrically short dipole located at the center of a cubic sample of WM as it is shown in Fig.|^ 
At / = 38 THz, the radiation intensity becomes enhanced 570 times compared to the absence of the WM. This 
enhancement was found as the radiation resistance of the dipole normalized to that calculated in free space. To get 
rid of dimensional resonances, we performed simulations for three different sizes of the sample: 12x 12, 16x 16 
and 20x20 nano wires increasing similarly the length of nano wires. The blue dashed curve corresponds to the 
analytically calculated enhancement of radiation for a transverse dipole in an infinite WM El. At the frequency of 
the topological transition, this value is 3.8 times lower than that for the axial dipole embedded in cubic sample of 
LiTa03 WM. We numerically checked - if our nanowires are PEC, the radiation enhancement vanishes for both 
axial and transversal dipoles. We also have done simulation for different values of the array period, wire radius, 
and different positions of the dipole in the cubic sample, which all confirmed our expectations. 

III. Conclusion 

In this work, we have pointed a topological transition in a medium of parallel nanowires from epsilon-negative 
material. Applying the spatially dispersive model we showed that the axial component of the effective permittivity 
considered as a function of both frequency and axial component of the wave vector vanishes at one value of this 
dual argument. This transition grants both elliptic and hyperbolic dispersion and in spite of rather high optical 
losses results in the significant (2-3 orders of magnitude) increase of the radiated power for a subwavelength 
dipole submerged in such the medium. It is important that this huge Purcell factor was achieved for an axially 
oriented dipole. We have done extensive full-wave numerical simulations with finite samples of the WM which 
confirmed our theoretical expectations. Notice, that we also studied the radiation of the submerged dipole into free 
space (radiative Purcell factor) which is also enhanced by our WM (by one order of magnitude) and demonstrates 
a specific pattern. These results will be shown in our presentation. 
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